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ABSTRACT 

A technique is presented for the numerical evaluation of the integrals 
occurring in the Rydberg-Klein-Rees method of calculating potential-energy 
curves. /fr 

INTRODUCTION 



The potential-energy curves for the bound states of diatomic molecules 

can be obtained from spectroscopic constants by using the Rydberg-Klein 

Rees method. Klein^ expressed the turning points of motion r and 

max 

r min -*- n terms of two auxiliary functions f and g: • 

r max = ^ + f /s) 1/2 + f 
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(f 2 + f/g) 1 / 2 - f 


Both of the functions f and g depend parametrically on the potential 
energy U and an additional parameter K = J(J + l)h 2 /2(j. where J is 
the rotational quantum number and p. is the reduced mass. The functions 
f and g are defined as 
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Here E(l, K) represents the vibrational-rotational energy levels, and 
I equals h(v + l/2) with ▼ representing the vibrational quantum number. 
From spectroscopic data the energy levels are generally expressible in 
the form 


P q. 

E(I,K) - V Y lm (r + 1/2) Z [J(J + l)] m (3) 

2=0 m=0 

where Y^ m are the usual spectroscopic constants. The upper limits of 
integration in (l) and (2) are obtained from E(l*,K) - U = 0. 

Rees° evaluated the integrals in (l) and (2) analytically for the 
special cases where E(l, K) was either quadratic or cubic in I, The 
results in the cubic case were not in a very convenient form for comput- 
ation. The results in the quadratic case, however, have been used as the 

4 5 

basis for rather extensive calculations by Vanderslice and coworkers ’ 

7 

and others ' . Molecules, whose data could not he adequately represented 
over the entire range hy a quadratic, were treated by piecewise fitting of 
quadratics to the energy levels. As pointed out hy Weisanan, Vanderslice, 

o 

and Battino this piecewise fitting can lead to errors. 

The difficulty in the numerical integration of (l) and (2) is caused 
by the fact that the denominator of the integrands has a zero at the upper 
limit of integration. JarmaJn^ cicumvented bhis problem in an^approxi- 
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mate manner by fitting [U - E(l, K)] 1 / 2 to an expression of the form 
c(l - l)“^. The constants c and a were evaluated by using two points 
very close to the upper limit of integration. Similarly, Se/Sk -was 
fitted to a quadratic in (I* - I). These approximations were then ana- 
lytically integrated to evaluate the contributions to f and g from 
regions close to I* . More recently, Weissman, Yanderslice, and 
Battino 8 introduced a new integration variable x = [U - E(l, K) ]^/ 2 in 
order to eliminate the singularity from the integrands. Although cor- 
rect, thiB procedure produces some unnecessary numerical inconvenience 
since the integrands are available as expressions in I and not the 
integration variable x. Thus, in numerical evaluation of the integrals 
one cannot use arbitrarily selected increments of x but rather must 
select increments of I and calculate increments of x as 

Ax = [U - E(I + AI, K)] 1 / 2 - [U - E(l, K)] 1 / 2 
The technique that will be described here does not require the fit- 
ting of the integrand near the upper limit, and further, it essentially 
retains the original integration variable. It is based on the fact that 
the singularities in the integrands of f and g can be easily removed 
by an integration by parts. 

ALTERNATE METHOD FOR NUMERICAL EVALUATION 
The quantity U - E that appears in the integrands of f and g is 
a polynomial of order p in the variable I with coefficients that de- 
pend on K. Introducing the notations x * v + l/2 and K = j(j + l), 
we have in dimensionless form 
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P 

(U - E)/E q = P(xj/c) = ^ A [( k ) x 1 

1=0 

where E^ is a constant with the dimensions of energy and where the 
coefficients in the polynomial P have the explicit form 

<1 

A i (,c) • (U6 10 • hm <m) ' E 0 

m=0 

In this notation the upper limit of integration I 1 ' in (l) and (2) 
corresponds to hx^(k) where x-j_(k) is the smallest, real, positive 
root of P(x;k). 

If the existence of the improper integral f is assumed, then (l) 
implies that the polynomial P(xjk) has a zero of order one at Xq_(k). 

If Kq is a particular value of k, then we can make a change in scale 
z = x/x^(Kq), and the polynomial P can "be written as 
P 

P(zj k ) ' A z (<)z Z = (1 - z)R(z;/c) - [1 - x 1 (k)/x 1 (/c 0 )]R(z;/c) (4) 

1^0 

where Aj(k) = x^(k^)A^(k) and where R(z ;k) is a polynomial in z of 
order p - 1. If K is chosen equal to then (4) takes on the 

simple form 

P 

P(z ;k) = A^(k)z 1 = (1 - z)R(zjk) (5) 

1=0 

where the subscript has been dropped from Kq for convenience, and where 

P-1 

l ' 'T' 1 

now Aj(k) = x 1 (/c)Aj (k) and R = / , B^(k)z . Comparing coefficients of 

f=0 

like powers of z in (5) establishes the relationship between and 
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A I as 


Xv*> 


Z “ 0, 1, • • • , p *■ 2 


B 2 (#c) =< 
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-A (k) 
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Z = p - 1 
P 


Also from (5) for z = 1 one can obtain Aj(k) = °* Combining this 

with (6) finally gives the result 


Z=0 


B z (k) =2 Aj(/c) Z = 0,1, . . ., p - 1 (7) 

J = Z+1 

The numerator of the integrand of g is essentially dp/d/c. This 
partial derivative is most conveniently evaluated by using the first 
equality in (4). After the differentiation is performed, K is again 
chosen to be equal to Kq» This gives, after dropping the subscript on 

V 

P 

dp/d/c s Q(z }K) = ^ C z (k)z 1 (8) 

1=0 


where 


Cj(k) = x^(k)c)Aj/c)k 

The integrals for f and g can be written in dimensionless form 
by using the Bohr radius a Q as a unit of length and Eq = fi 2 /2ua.Q. 
Also, making a change in scale to the new variable z and considering 
the particular value k that is equal to Kq result in 
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where the subscript on k q has again been deleted for convenience. 
Substituting (5) for P(z;/c) in (9) and (10) and performing an integra- 
tion by parts give 



In (11) and (12) we have achieved the objective of removing the singu- 
larity of the integrands at the upper limit. For the important case 
< = 0 , the polynomial coefficients required to evaluate the integrands 
of (ll) and (12) have the relatively simple forms 


P 



No difficulties were encountered in the numerical application of (ll), 
(12), and (13) to the calculation of potential- energy curves for several 
diatomic molecules by using standard integration techniques. In fact. 



7 


the results of single and double precision calculations for the ground 
state of Hg agreed to at least seven figures when using the spectro- 
scopic constants of Weissman, Vanderslice, and Battino®. 
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